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ABSTRACT 

We present an analysis of the early Kepler observations of the previously discovered transiting planet 
HAT-P-7b. The light curve shows the transit of the star, the occultation of the planet, and the orbit 
phase-dependent light from the planet. In addition, phase-dependent light from the star is present, 
known as "ellipsoidal variations" . The very nearby planet (only 4 stellar radii away) gravitationally 
distorts the star and results in a flux modulation twice per orbit. The ellipsoidal variations can confuse 
interpretation of the planetary phase curve if not sclf-consistently included in the modeling. We fit the 
light curve using the Roche potential approximation and derive improved planet and orbit parameters. 
Subject headings: binaries: eclipsing 



1. INTRODUCTION 

Kepler is a reconnaissance mission to obtain time-series 
optical photometry of ~150,000 stars in order to de- 
termine characteristics of Earth-size and larger extra- 
solar planets: frequencies, sizes, orbital distributions, 
and correlations wit h the properties of the host stars 
(jBorucki et al.|[2010( ). To achieve these goals, Kepler re- 
quires exceptional lon g-term photometric stability and 
precision (Koch et al. 20101) . In addition to the discov- 
ery aspect of the mission, this unprecedented photomet- 
ric capability provides exquisite observations of a host 
of astrophysical objects, including the previously known 
extrasolar planets in the Kepler field of view. In this 
Letter, we examine the early Kepler observations of the 
planet HAT-P-7b. 

Discovered via the HATNet project, the transiting 
planet HAT-P-7b revolves in a tight circular orbit 
(a=0.038 AU, P=2. 20473 d) a round a bright (V=10.5 
mag) F6 star (iPal et all 120081). T he proximity to its 
6350 K host star ( Pal et aTT 20081 ) means the planet 
is highly irradiated, resulting in very high tempera- 
tures (^2140 K), m aking it an extreme pM-type planet 
(|Fortnev et al.|[200S ). Using observatio ns of t he Rossiter- 
McLaughlin eff'ect. iNarita et al.l ()2009f ) and iWinn et all 
12009) find that the planet's orbital axis is extremely 
tilted compared to the star's spin axis, and probably even 
retrograde, implying an interesting formation and orbital 
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evolution history. Additionally, the radial velocities ex- 
hibit an acceleration, suggesting the presence of another 
body in the syste m, perhaps responsible for the tilted 
orbit (jWinn et al.l l2009V 

Observations of HAT-P-7 during the 10 days of 
commissioning of the Kepler photometer revealed the 
presence of an occultation (also known as a "sec- 
ondary eclipse") as the planet passes behind the star 
(jBorucki et al.l 120091 ). These data also show the phase 
"refiected" light from the planet, a result of both scat- 
tered and thermal emission. Additional observations of 
HAT-P-7 have been obtained, and we present these in §2. 
In §3 we describe our modeling method, with particular 
emphasis on the ellipsoidal variations from the star. In 
§4 we present and discuss our findings. 

2. OBSERVATIONS 

HAT-P-7 was monitored continuously for 33.5 d dur- 
ing the 2009 May 13-Jun 15 "Quarter 1" (Ql) epoch 
in short-cadence mode. The 42 Kepler CCDs were read 
out every 6 s and co-added on-board to achieve approx- 
imately 1 min sampling cadence. The photometer has 
no shutter, so an overscan region is used to remove the 
effects of smearing during readout. In 6 s exposures 
HAT-P-7 saturates the CCD; however, because the Ke- 
pler photometer is such a stable platform, this does not 
hamper the relative precision and superb photometry 
is possible. For details on the design and perform ance 
of the Kepler phot ometer, see iKoch et al.l (|2010[ ) and 
iJenkins et ahl (pOlOl) . 
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Fifteen complete transits of HAT-P-7 were observed, 
and after removing 21 cosmic rays via a +4 sigma rejec- 
tion from a 30-min wide running median, 49,023 mea- 
surements remained. The background-subtracted light 
curve exhibited a gradual drop in counts (0.13%), sus- 
pected to be due to a drift in focus coupled with a fixed 
extraction aperture. This trend was removed using a 
cubic polynomial (masking out the transits). Because 
we are interested in phenomena on the orbital timescale 
(2.2 d), we chose this modest detrending to leave in as 
much power as possible on orbital timescales. Fluxes 
were normalized at phase 0.5 (mid occultation) because 
the planet is hidden at this phase and we see starlight 
alone. Sho rt-timescale systema tic calibration features 
are present ([Gilliland et al.|[20Tol) . but they does not af- 
fect the analysis other than being a noise term. Uncer- 
tainties were estimated by taking the median of the rms 
deviations in 30-min bins, resulting in 150 ppm per 1-min 
datum. 

The detrended and phase-folded light curve is shown 
in Fig. 1, where the upper panel includes our fit to the 
transit, and the lower panel highlights the eclipse of the 
planet and the phase-dependent light from the planet 
and star. Notice that the maximum does not occur just 
outside occultation as one would expect for simple "re- 
flection" from the planet. Two maxima occur 0.15-0.20 
away in phase, a result of the "ellipsoidal variation" of 
the star's light, as discussed in §4. 

The absolute timing is still preliminary in this early 
version of the data calibration pipeline, but relative tim- 
ing precision is reliable (though without the modest B JD 
correction, which is unimportant for our investigation 
over the 33.5 d time series). So we do not report the 
value for the epoch of transit TO, though of course it is 
a parameter in the model fitting. Also, the period mea- 
sured is based only on these 33.5 d, so the precision is 
not as high as would be if other epochs many cycles away 
were included. 

3. MODELING 



We employ the ELC code of Oro sz &: Hauschildtl 
(|2000f ) to model the transit, occultation, and phase- 
varying light from the planet and star. The code si- 
multaneously fits the photometry along with several ob- 
servational parameters: the radial velocity K, and the 
mass and r adius of the star. The K amplitude was 
taken from I Winn et al.l (|2009l ). and the mass and ra- 
dius from the asteroseismology analysi s of the Kepler 
data (jChristensen-Dalsgaard et 311120101) . These param- 
eters are not fixed; rather the models are started and 
steered toward them via a chi-square penalty for devia- 
tions. Markov chain Monte Carlo and genetic algorithms 
were used to search parameter space, find the global chi- 
square minimum, and determine confidence intervals of 
the fitted parameters. 

The analytic model of iGimenej ()2006fl for a spheri- 
cal star and planet is not sufficient to model the phase 
variations. Therefore ELC is used in its full numeri- 
cal mode: the star and planet are tiled in a fine grid, 
and the intensity and velocity from each tile is summed 
to give the light curve and radial velocities. Limb and 
gravity darkening are included, and the gravitational dis- 
tortions are modeled assuming a standard Roche poten- 
tial. We employed a blackbody approximation evalu- 



ated at a wavelength of 6000 A, and a hybrid methocQ 
where model atmospheres are used to determine the in- 
tensities at the normal for each tile and a parameterized 
limb darkening law is used for other angles; the model is 
then filtered through the Kepler spectral response func- 
tion (spanning roughly 4250-8950 A, peaking at 5890 A 
with a me an wavelength of 6400 A — see iKoch et al.l 
IpoToh and lVan Cleve & Caldwell (|2009t) . The two meth- 
ods give essentially the same results, with the exception 
of a higher albedo from the hybrid models. Interest- 
ingly, the blackbody models yielded significantly lower 
chi-square values, so we quote the blackbody model val- 
ues in this work. This needs to be kept in mind when in- 
terpreting the temperatures and albedos estimates given 
in §4. We use a 2-parameter logarithmic limb darkening 
law and adopt an eccentricity of zero consistent with the 
radial velocities and phase of occultation. We assume 
the planet is tidally locked in synchronous rotation. For 
more details on us i ng EL C to model exoplanet data see 
iWittenmv er et al.l (|2005[ ). _ 

Following the prescription of iWilsonI ()1990D . the light 
from the planet is modeled as the sum of an isothermal 
component (with temperature Tp that essentially adds a 
constant flux at all phases outside of occultation), and 
a "reflection" component on the day hemisphere. The 
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local temperature is given by 

which comes from assuming the fluxes and temperature 
are coupled by the Stefan-Boltzmann law. The bolomet- 
ric albedo Aboi, also known as the "heat albedo" is the 
ratio of re-radiated to incident energy, and should not 
be confused with the Bond albedo. For stars, a radiative 
atmosphere has Aboi — l^ (local energy conservation) while 
for convective atmospheres a value of 0.5 is appropriate 
(h alf the energy is rad iated, half gets redistributed) — 
see lKallrath fc Milond ([2999) for a full description of the 
method. In our model we allow Ai,oi and the tempera- 
ture ratio Tn/T^, to be free parameters, keeping fixed 
at 6350 K (jPal et al.1120081 ). The term in brackets is the 
local refiection factor TZ, equal to the ratio of the total 
radiated flux (internal -|- re-radiated) to internal flux. 
Thus TZ > I, with 7?, = 1 on the night side. For an iso - 
lated planet, Tp is very low (e.g. iBurrows et al.l ()2006f ) 
assume 50 K), but for an irradiated planet, much of the 
incident energy is eventually re-distributed, bringing the 
night-time temperature up to much higher temperatures. 

In the ELC model, the local emitted flux is completely 
dependent on the local temperature. The local tempera- 
ture depends on the mean effective temperature, the local 
gravity, and irradiation. The irradiation term is the sum 
over all the visible tiles on the star as seen from each tile 
on the planet, and includes the ellipsoidal shape of the 
star, gravity darkening, limb darkening, and penumbra 
correction (accounting for the fact that parts of the star 
are not visible because they are blocked by the local hori- 
zon). Note that given the close proximity of the planet 
to the star, no part of the planet sees a full hemisphere 
of the star. There is no explicit scattering term in the 

^ ELC can also fully employ stellar atmosphere intensities, 
where no parameterized limb darkening is used. But preliminary 
tests gave significantly worse fits unless the stellar temperature was 
allowed to be several thousand degrees hotter. The cause seems to 
be related to the very wide Kepler bandpass and the lack of free- 
dom to adjust the limb darkening. 
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reflection: the incident radiation heats up the planet and 
is always re-emitted locally. Scattered light is implicitly 
accounted for by allowing the global temperature Tp to 
be de-coupled from the day-side temperature. 

The ELC model provides a good fit to the observations. 
The chi square is 57,225 for 49,016 degrees of freedom 
(reduced x^=:1.167). Uncertainties on the parameters 
were boosted by a factor of a/xB" (=8%) to account for 
the formally high x^, which we attribute to systematic 
non-Gaussian noise. The values of the parameters are 
listed in Table 1. 

4. RESULTS AND DISCUSSION 

The 1.8 Mjup planet orbiting only 4.1 stellar radii from 
its host star induces a tidal distortion on the star, chang- 
ing its shape from oblate (not spherical, due to its rota- 
tion) to a more triaxial shape, with the longest axis along 
the direction toward the planet and the shortest axis per- 
pendicular to the orbital planeO This shape causes the 
well-known "ellipsoidal variation" effect seen in binary 
stars: a modulation in light at half the orbital period, 
with maxima at phases near 0.25 and 0.75, and unequal 
minima at phases 0.0 and 0.5. The ellipsoidal variation 
is primarily a geometrical (projected surface area) effect 
whose relative amplitude depends on the mass ratio and 
the inclination of the binary system. Given the tight con- 
straint on the inclination from the eclipses, they may pro- 
vide some limits on the mass ratio, and hence the mass 
of the planet. In the optical, where the phase-dependent 
planet-to-star flux ratio is small, the presence of the stel- 
lar ellipsoidal variation becomes significant. Neglecting 
its contribution can lead to a confused interpretation of 
the phase curve and thus an incorrect measurement of the 
albedo and phase-dependent scattered/thermal emission. 

The presence of ellipsoidal variations in e xoplanetary 
system s was antic ipated bv iLoeb fc Gaudil (2003) and 
[Drake (2003'), and iPfahLsLalJ (12008') present a detailed 
theoretical investigation of the tidal force on the star by 
the planet. Figs. 1 and 3 show the first detection of 
the effect in an exoplanet system. In Fig. 3 we show 
the light curve binned to 30 min, along with the best- 
fit model decomposition. The dotted curve shows the 
stellar-only ellipsoidal light curve while the dashed curve 
is the planet-only light curve (offset vertically) . The solid 
curve is the sum of the two and equals our best-fit model. 
The planet-only model can match the light curve only 
very near phase 0.5; it is a very poor fit to the data at 
other phases, indicating the need for the ellipsoidal vari- 
ation component. The amplitude of the ellipsoidal com- 
ponent is 37.3 ppm, detectable only because of Kepler's 
high-precision photometric capability. 

Ellipsoidal variations arise as a consequence of gravity 
on a luminous fluid body. Within the Roche framework, 
its amplitude is exactly known if the mass ratio, incli- 
nation, and stellar radius are known. However, since 
in exoplanet systems the star is not expected to be in 
synchronous rotation, the Roche potential is an approx- 
imation, albeit a good one since the maximum tidal dis- 
tortion AR/R is only 10""*. However, we stress that the 
model presented here is only a starting point, based on 

^ More exactly, in a Roche potential there are 4 distinct radii: 
towards the companion, perpendicular to the orbit plane, along the 
direction of motion, and away from the companion. 



the well-developed and suc cessful Roche mod el. Effects 
such as those discussed in iPfahl et all ()2008l ) based on 
the equilibrium tide approximation can be present. As 
more Kepler data become available it will be interest- 
ing to try to distinguish between these approximations, 
as they can lead to measuring internal properties of the 
star. For example, in the frame of the rotating binary, 
the star is spinning retrograde if Pspin < Port, which is 
expected to be true for most hot Jupiter planets. This 
could induce a phase shift in the ellipsoidal variations, as 
the star's spin "drags" the tide away from the line of cen- 
ters. If a phase-lag not associated with the planet light 
can be unambiguously measured, it may be possible to 
put constraints on the tidal Q factor. This asymmetry 
could also lead to interesting orbital dynamics effects. 
The detailed shape and amplitude of the ellipsoidal vari- 
ation also depends on the stellar envelope (convective vs. 
radiative - Pfahl et al. (2008)), thus potentially offering 
a probe of stellar interiors. However, for HAT-P-7, the 
situation is complicated by the f act that the spin axis 
is not aligned with the orbit axis mn eFaH [20091) . so 
these optimistic statements must be tempered with cau- 
tion. 

Fig. 2 shows model light curves that illustrate the con- 
tributions of the star to the optical light. The lower curve 
is the stellar light only, showing the ellipsoidal variations. 
It is not constant outside of transit, and its amplitude 
can be significant compared to the depth of the occulta- 
tion. Above this are four curves that include the planet 
contribution. The lowest is for an isothermal isotropic 
planet with no phase-dependent scattered or re-radiated 
emission. To first order, it is essentially an additive offset 
to the star's light, the amount depending on the relative 
radius and temperature of the planet. If the absorbed in- 
cident energy from the star is not perfectly uniformly re- 
distributed around the planet, the day side of the planet 
will be hotter than the night side, resulting in an orbit- 
phase dependent modulation peaking at the sub-stellar 
point (phase 0.5). Or, if the atmosphere scatters the 
incident radiation, again the day side will be brighter 
than the night side (though not necessarily hotter) . The 
upper three curves show this phase-dependent effect for 
different values of the heat albedo, which is a proxy for 
scattered and re-radiated emission on the day side. Note 
that advection can move the peak downwind and pro- 
duce a phase shift in the planet's emission, as seen in the 
infrared phase curve of HD189733 (Knutson et al. 2007); 
this is not included in these models. 

Given a stellar temperature estimate, the ELC model 
can yield day and night hemisphere temperatures. Be- 
cause the planet is not black it contributes light at all 
phases (other than occultation) and its flux contribution 
at each orbital phase, relative to ellipsoidal and other ef- 
fects, allows its temperature to be estimated. In partic- 
ular, the ability to measure relative temperatures arises 
from the requirement of getting the occultation and tran- 
sit depths correct given the tight geometric constraints 
(e.g. a 2600 K planet reduces the 6000 A transit depth by 
29 ppm compared to a zero temperature planet). We find 
the peak temperature on the planet (at the sub-stellar 
point) is 3160 K, but a more meaningful flux-weighted 
day temperature is 2885 ±100 K. This is notab ly higher 
than the 2560 ±100 K estimate bv lBorucki eTalJ ([20091). 
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We caution that the day-side temperature estimate is de- 
rived from the assumption that the light is entirely ther- 
mal in origin with no scattered component; thus this is 
a maximum day-side temperature estimate for the given 
heat albedo. We find the night side temperature Tp is a 
surprisingly hot 2570 ±95 K, only 590 K less than the 
peak day temperature. 

As seen in Fig. 3, the depth of the occultation is 
85.8 ppm, and the night-side contribution of the planet 
is 22.1 ppm, giving a day-side peak flux enhancement of 
63.7 ppm. Combining these flux ratios, temperatures, 
and the relative radii, we can atte mpt to estini a te the 
geometric and Bond albedos — see iRowe et al.l pOOl 
for a discussion. The planet's albedo is key to the en- 
ergy balance in the planetary atmosphere and its im- 
portance for the d ayside emissi o n is d is cussed in sev- 
eral p a pers (e .g. Seagcr ct al. (200^, IBurrows et al.l 
(|2006D. iLopez^Morales fc Seager ( 20071) ). The geomet- 
ric albedo is given by the ratio of planet to star flux 
at phase 0.5, scaled by the ratio of surface areas: 
Ag = {Fp / F^) / [Rp / . Using the peak planet to star 
flux of 63.7 ppm yields Ag=0.18. This agrees with 
the <0.20 value found fo r the pM-clas s planet CoRoT- 
Ib (|Snellen et al.l (|2009f ). lAlonso et al.l (l2009b)), but is 
significantly higher than s ome other p lanets, notably 
Ag <0.Q8 for HD 209 458b (IRowe et al.([^08fl . and 0.06 
±0.06 for CoRoT-2b ()Alonso et al.ll2009aD . The equilib- 
rium temperature and Bond albedo are related by the 
definition: = Tt{R*/2af[f{l - Ab)] where / is the 
redistribution factor and equals 1 for complete redistri- 
bution or 2 if the incident flux is re-radiated only on the 
day hemisphere. Setting /=1 and equating the night-side 
temperature with the equilibrium temperature should in 
principle allow one to solve for the Bond albedo since all 
the other terms are known. However, this fails because 
the maximum equilibrium temperature possible (when 
Ab—Q) is 2213 K, considerably less than the ELC esti- 
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mate. It may be that our 6000 A brightness tempera- 
ture estimate exceeds the equilibrium temperature sim- 
ply because the planet is not a blackbody and at other 
wavelengths lower temperatures may be measured. We 
speculate on two other possibilities, both related to the 
presence of a third body in the HAT-P-7 system. First, 
the planet may genuinely be hotter then the equilibrium 
temperature due to non-radiative heating, perhaps tidal 
heating due to an encounter with another object (recall 
the >8 6 degree offs et between the orbit and stellar spin 
axes — iWinn etaLl (2009), Narita ct al. (2009)). A more 
mundane, but perhaps more likely, explanation is that 
light from a third body is contaminating the flux ratios. 
The acceleration term in the radial velocities suggests 
the presence of another body (jWinn eFan [20091) . When 
more Kepler transits are available, we can check for tran- 
sit timing variations and address the issue of potential 
third light contamination. 

In closing, the Kepler light of HAT-P-7 curve reveals 
ellipsoidal variations with an amplitude of approximately 
37 ppm. This is the first detection of ellipsoidal varia- 
tions in an exoplanet host star, and shows the precision 
Kepler is capable of producing at this early stage. For 
comparison, a transit of an Earth-analog planet around 
a Sun-like star would produce a signal depth of 84 ppm, 
a factor of 2 larger than this effect. 
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Fig. 1. — upper: Detrended and phase-folded light curve at 1-min cadence along with the ELC fit. For scale, the horizontal bars show the 
vertical size of the lower panel, lower: Light curve averaged in 5 min and 75 min bins. The double-humped shape is due to the ellipsoidal 
variations of the star plus light from the planet. 
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Fig. 2. — ELC model light curves for half an orbit. The bottom curve is the stellar light curve exhibiting the ellipsoidal variations. Above 
this are planet+star light curves with increasing bolometric (heat) albedos: 0.0, 0.3, 0.6 and 0.8. If the albedo is 0.0 the planet emits a 
constant amount of light resulting in a simple offset in flux. As the albedo increases, the light from the planet becomes more pronounced. 
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Fig. 3. — Phase-folded Kepler light curve with best fit model (solid curve). Also shown are the component stellar-only ellipsoidal model 
(dotted) and the planet-only model (offset by +1; dashed). Both the data and models have been cast into 30-min bins. 
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TABLE 1 
HAT-P-7 System Parameters 



Psyrctmctcr 


Value 


Uncertainty 


Unit 


fji a 


6350 


— 


K 




1.53 


0.04 


Mq 




1.98 


0.02 


Rr^ 




212 


5 


m 


Orbital inclination, i 


83.1 


0.5 


degrees 


Orbital period, P 


2.204733 


0.000010 


days 


Star-to-planet radius, R»/Rp 


12.85 


0.05 




limb dark coefficient x 


0.58 


0.08 




limb dark coefficient y 


0.21 


0.13 




Mass of planet, Mp 


1.82 


0.03 


Mjup 


Radius of planet, Rp 


1.50 


0.02 


Rjup 


Semimajor axis, a 


8.22 


0.02 


R© 


Bolometric (heat) albedo, Ajoj 


0.57 


0.05 




Tp (night side) 


2570 


95 


K 


Tp (average day side) 


2885 


100 


K 



Note. — 

Held fixed at 6350 K: IPal e t aT (WOS) 
^ Steered toward 1.52 ±0.036: Christcnse n-Dalsgaard et all 1 12010) 
^ Steered toward 1.991 ±0.018; Cliristensen- Dalsgaard et all 1 12010ft 

Steered toward 211.8 ±2.6; iWinn et al.. t200"9h 



